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We discuss a number of basic physical mechanisms relevant to the formation of the proximity 
effect in superconductor/normal metal (SN) systems. The physics of this effect, most notably the 
phenomenon of density of states depletion in normal metals adjacent to superconductors, sensi- 
tively depends on various aspects of quantum interference and chaotic charge carrier dynamics in 
mesoscopic systems. Specifically, we review why the proximity effect sharply discriminates between 
systems with integrable and chaotic dynamics, respectively, and how this feature can be incorpo- 
rated into theories of SN systems. Turning to less well investigated terrain, we discuss the impact 
of quantum diffractive scattering on the structure of the density of states in the normal region. We 
consider ballistic systems weakly disordered by pointlike impurities as a test case and demonstrate 
that diffractive processes akin to normal metal weak localization lead to the formation of a hard 
spectral gap - a hallmark of SN systems with chaotic dynamics. Turning to the more difficult case 
of clean systems with chaotic boundary scattering, we argue that semiclassical approaches, based on 
classifications in terms of classical trajectories, cannot explain the gap phenomenon. Employing an 
alternative formalism based on elements of quasiclassics and the ballistic cr-model, we demonstrate 
that the inverse of the so-called Ehrenfest time is the relevant energy scale in this context. We discuss 
some fundamental difficulties related to the formulation of low energy theories of mesoscopic chaotic 
systems in general and how they prevent us from analysing the gap structure in a rigorous manner. 
Given these difficulties, we argue that the proximity effect represents a basic and challenging test 
phenomenon for theories of quantum chaotic systems. 
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I. INTRODUCTION 

Superconductors attached to normal metals tend to 
export aspects of their anomalous properties into the ad- 
jacent normal metal region. Termed the proximity effect, 
this tendency manifests itself in phenomena such as the 
Josephson effects, anomalous charge and thermal trans- 
port coefficients, and altered thermodynamic properties 
(see e.g. Refs. |l]-|3| for extensive reviews). All these 
effects find their common origin in the fact that a super- 
conductor pairing field amplitude may penetrate deeply 
into the normal metal region, before decaying due to the 
absence of attractive interaction mechanisms. 

Perhaps the most direct manifestation of the presence 
of a finite pairing field amplitude in a normal metal adja- 
cent to a superconductor is a massive suppression of the 
single particle density of states (DoS). This phenomenon, 
first noticed almost forty years ago, has a long history of 
research. Beginning with a classic paper of de Gcnncs 
and Saint- James g, the normal metal DoS of a wide 
family of superconductor-normal-metal hybrid systems 
(SN-systems) has been investigated. Among these are su- 
perconductors coated with normal metal thin films, both 
clean and disordered bulk diffusive compounds in 

contact with superconductor terminals [ |_ lOfl a nd various 
types of SN quantum dot structures 11 ~[h| . 

In general the quantitative functional form of the en- 
ergy dependent DoS, z/(e), of any of these systems de- 
pends on the system category, and on a number of rel- 
evant material parameters. Interestingly, however, a 



simple ordering principle exists, as first pointed out by 
Melsen et al. p^ , |l3] |, according to which the qualitative 
structure of the DoS falls into one of only two differ- 
ent categories: For systems whose normal compound is 
integrable - in the sense that the dynamics of the sin- 
gle particle degrees of freedom becomes integrable in the 
classical limit - the DoS vanishes in a close to linear fash- 
ion P, pli||l3| ] as the excitation energy e approaches zero. 
Of course, due to the inevitable presence of imperfections, 
truly integrable systems in condensed matter physics do 
not exist. However, a number of systems, e.g. thin clean 
metallic films, come reasonably close to the integrable 
limit implying that the proximity DoS obeys a continu- 
ously decreasing behavior down to very small excitation 
energies e. In contrast, the DoS of genuinely chaotic sys- 
tems - significantly disordered systems or chaotic quan- 
tum dot structures - is gapped: v(e) falls to zero at some 
finite threshold energy e*. For energy values |e| smaller 
than e*, the DoS is strictly vanishing. 

Beyond its practical relevance in the physics of SN sys- 
tems, this twofold categorization scheme is of noticeable 
conceptual interest. Indeed, the observation that a quan- 
tity as elementary as the single particle DoS of may serve 
as an indicator sharply discriminating between integrable 
and chaotic dynamics has led to a wave of interest in the 
interplay between superconductor correlations and quan- 
tum chaotic dynamics |ll| [D^ , |l7| -pl| . Unfortunately, the 
established tools of theoretical mesoscopic superconduc- 
tivity - quasiclassical Green function methods and dia- 
grammatic techniques - are not straightforwardly appli- 



1 



cable to clean chaotic SN structures. As an alternative to 
direct numerical analyses, an approach based on classical 
path counting, augmented by a Bohr-Sommerfeld quan- 
tization rule, has become popular. Referring to a more 
substantial discussion to below, we here merely note that 
the Bohr-Sommcrfeld (BS) approach readily produces a 
massive suppression of the low energy DoS of chaotic sys- 
tems fl7],|l8| (as compared to the integrable case.) Fur- 
thermore, the energetic extension of the depletion region 
turns out to be in good accord with the results obtained 
by direct diagonalization p^]l3] , |l^ , |l8| , f2"l[ . In this sense, 
the BS quantization scheme represents a comparatively 
simple and widely applicable tool for exploring the DoS 
of SN systems. 

There is, however, one key feature of the DoS gap that 
is missed by the BS-scheme: for chaotic systems the gap 
tends to be hard (in the sense that v(e) strictly vanishes 
for energies |e| < e*.) Hard gaps are found for a wide class 
of chaotic SN systems: For the particular case of diffusive 
SN systems, the structure of a hard gap has been ob- 
tained analytically by quasiclassical methods p2],|6|-^0|. 
A numerical study by Pilgram et al. ]23[ | has shown that 
the gap of weakly disordered systems, i.e. systems whose 
spatial extent is smaller than the elastic mean free path, 
is also hard (within, of course, the numerical accuracy). 
The same applies to the numerically obtained DoS pro- 
files of systems with smooth chaotic boundary scattering 
displayed in Ref. jl7|l . Finally, chaotic systems that are 
completely ergodic, in the sense that their quantum me- 
chanics can be modelled by random matrix theory, also 
exhibit a hard gap |L2 13 1. 

In contrast, the very construction of the BS scheme im- 
plies that it categorically produces a finite (albeit small) 
DoS at non-zero energies |l|,|l|,0,|8U|^] . That no gap 
is seen has to do with the fact that the the BS approach 
does not fully account for the quantum mechanical na- 
ture of the quasi-particle propagation in an SN environ- 
ment. More specifically, the elimination of all states be- 
low a certain threshold energy is caused by mechanisms 
of quantum diffraction similar to those responsible for 
weak localization corrections in normal metals. The es- 
sential difference to the weak localization processes op- 
erative in normal mesoscopic systems is that diffraction 
processes in SN systems are not small. Indeed it is one of 
the main objectives of the present paper to show that at 
low enough energies these 'corrections' diverge and that a 
non-perturbative treatment is needed to correctly obtain 
the gapped DoS. 

Qualitatively, the results of our analysis may be sum- 
marized as follows: in chaotic SN systems, there are at 
least two different time scales that are of relevance for 
the low energy DoS. The first is the average classical 
time t c it takes to establish contact with the supercon- 
ductor. This scale can be set by be the diffusion time 
time through an extended disordered system, the inverse 
tunneling rate through a nearly insulating superconduc- 
tor normal metal interface, or the average ballistic time 
of flight through a clean quantum dot. The second time 



scale is, in general, of quantum mechanical nature, and 
is the time id a minimal wave package takes to lose mem- 
ory of its classical initial conditions. In a medium with 
hard (s-wave) scatterers, will coincide with the classi- 
cal scattering time r. In systems with smooth scattering 
potentials - no matter whether disordered or clean - 
is determined by the Ehrenfest time £e [0, that is to 
be discussed extensively below. For the structure of the 
low energy DoS, the larger of these two timescales is of 
foremost relevance. For example, in a typical diffusive 
system with hard impurities, the scattering time is much 
shorter than both the diffusion time and the tunneling 
time into the superconductor. In this case, t c > and 
the parametric scale determining the gap edge p|,p|^10[ is 

L c ■ 

In this paper the focus will be on the complementary 
class of systems with td 3> t c . We will begin by con- 
sidering the prototypical class of ballistic systems weakly 
polluted by some s-wave impurities and well coupled to a 
superconductor. For such systems the scattering time is 
larger than the ballistic time of flight (and the tunneling 
time to the superconductor) and the above condition is 
met. Operating within a synthesis of the quasiclassical 
formalism and the ballistic non- linear cr-model [^5|,^6j , we 
will identify the relevant quantum diffraction processes, 
and show how these contributions lead to a vanishing of 
the DoS for energies e < e* ~ h/(2r). 

Building on the insight gained in this analysis, we 
will then proceed to the more complex problem of clean 
chaotic systems with smooth potentials. In this case, 
the relevant corrections to the classical picture emerge 
from an interplay of quantum wave packet spreading and 
the exponential trajectory divergence inherent to chaotic 
dynamics. The characteristic time scale for these pro- 
cesses is the aforementioned Ehrenfest time. Below we 
will demonstrate that the Ehrenfest time is the relevant 
energy scale for the DoS of chaotic systems with smooth 
scattering. It should be mentioned right away that the 
quantitative parts of this part of the analysis will include 
certain assumptions and, therefore, do not attain the sta- 
tus of a rigorous calculation. Further, our formalism is 
not capable of exploring structures in the close vicinity of 
the estimated gap edge ~ h/t^. Notwithstanding these 
deficiencies, our analysis demonstrates the suppression 
of the DoS below its semiclassically estimated value. Al- 
though we cannot mathematically prove that this quan- 
tum suppression mechanism produces a hard gap, the 
close structural analogy with the hard scattering prob- 
lem makes us believe that this is the case. That a hard 
gap of the order of the inverse Ehrenfest time should be 
generated has already been conjectured by Lodder and 
Nazarov 0. 

We would like to point out that the phrase 'hard gap', 
as used in this paper, needs some qualification: refined 
analysis of the spectral structure of diffusive Andreev sys- 
tems has shown the existence of a tail below the quasi- 
classical gap edge Q that is exponentially small in g, 
where g ^> 1 is the dimensionless conductance of the 
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normal component of the system. Similar contributions 
have been indentified in billiards with randomly fluctuat- 
ing geometry In diffusive systems, the formation of 
these anomalous 'tail' states is due to exponentially rare 
potential configurations which lead to an effective decou- 
pling of the superconductor and normal metal jlo| . An 
interesting question, not addressed in the present paper, 
is whether similar states exist in clean chaotic environ- 
ments. It should be emphasized, though, that tail state 
formation is a mechanism conceptually independent from 
the quantum diffractive suppression of the semiclassical 
results discussed in this paper. 

Finally, let us mention that much of the conceptual 
ground on which our analysis of the DoS in chaotic sys- 
tems with smooth boundaries is based was laid in a sem- 
inal paper by Aleiner and Larkin J24|] (AL). In this work 
the fundamental importance of the Ehrenfest time for 
the formation of weak localization corrections in normal 
chaotic systems was realized and substantiated. Many of 
the ideas first formulated by AL have been included in 
our analysis of the superconductor problem. 

The paper is organized as follows. In section ||, we 
discuss some general physics associated with Andreev re- 
flection in SN structures, introduce the BS approach and 
discuss the role of quantum diffraction processes in estab- 
lishing the gap in the spectrum of the normal region. In 
section III we discuss how the BS approach may be recov- 



ered from the quasiclassical Eilcnbcrger equation in the 
diffractionlcss limit, and establish a simple trial solution 
of the quasiclassical equations in the presence of a col- 
lision operator. In section IV we introduce a functional 



integral formulation for the disordered Green function of 
an SN system in terms of a ballistic sigma-model. In 
section |\] we show how this formalism may be applied to 
determine the spectrum of a SN system with bulk, ballis- 
tic s-wave disorder, and establish the presence of a hard 

fof the order of the inverse scattering time. In section 
we turn to the study of the DoS in clean chaotic SN 
systems, in the presence of a smooth scattering potential, 
and discuss the appearance of the inverse Ehrenfest time 
in the structure of the DoS. Section VII concludes with 
a discussion and summary. 



II. QUALITATIVE PHYSICS OF THE GAP 

The basic microscopic mechanism of coupling between 
a superconductor and an adjacent normal metal is An- 
dreev reflection p7|] : a normal metal electron impingent 
on the superconductor interface is converted into a hole 
(charge balance being restored by the formation of a su- 
perconductor Cooper pair) which is reflected back into 
the normal metal. This type of scattering has a number 
of peculiar properties which are at the root of the forma- 
tion of the DoS gap and which are summarized below for 
later reference: 

• Andreev reflection is a retro-reflection, i.e. the hole 



is essentially reflected back into the path of the in- 
cident electron where 



• there is a slight angular mismatch set by the ex- 
citation energy e of the electron above the Fermi 
energy Ep. 



• For an electron of energy e, the reflected hole has 
energy — e. 



• the hole acquires a scattering phase, which takes 
the value of it/2 — ip (for an excitation energy e 
far below the superconducting gap) , where <p is the 
value of the superconducting order parameter at 
the point of incidence. 



From these properties, the origin of the DoS suppres- 
sion in normal metals adjacent to superconductors may 
be readily understood. Assume that an electron of some 
small excitation energy e is injected somewhere into the 
normal metal (at the beginning of the solid line in Fig. |l| 
say.) The electron will travel along some path 7 in the 
normal metal and at some point hit the superconductor 
where it is converted into a hole. The hole then propa- 
gates approximately (up to the above mentioned angular 
mismatch) back along the path 7 (dotted line in Fig. [I]) 
and eventually arrives in the vicinity of the point of in- 
jection. Seen as a whole, this process amounts to the 
formation of a local net pair amplitude ~ (i/^ (e)V^ (e)) , 
i.e. an amplitude for the creation of an electron and an 
annihilation of a hole, in the normal metal. Importantly, 
this induced Cooper pair amplitude is not weighted by 
rapidly fluctuating quantum mechanical phases. To un- 
derstand this point, notice that the quantum phase of 
the incident electron is essentially given by the classi- 
cal action 5* 7 (e) for traversal along 7. The action of the 
Andreev reflected hole (a missing electron) is given by 
— S 7 (— e). These two contributions largely cancel against 
each other, implying that a slowly fluctuating net ampli- 
tude obtains. 

Unitarity implies that the formation of a non- vanishing 
Cooper pair amplitude ~ (il>(e)'ip(e)) must be com- 
pensated for by a diminished 'normal state' amplitude 
(ip' (e)ip(e)) or, in other words, by a diminished local sin- 
gle particle density of states (see below for a more pre- 
cise phrasing of this argument.) It is also clear that this 
suppression mechanism can only be effective for small 
energies e, and for spatial regions not too distant from 
the superconductor: For increasing e or, equivalently, in- 
creasing length of the characteristic paths 7, the action 
difference <S 7 (e) — S' 7 (— e) increases, electron and hole run 
out of phase and the induced pair amplitude randomizes 
to zero. 
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pling to a superconductor (dark shaded region) leads to the 

formation of a stable pair amplitude ~ (^xj}) in the normal FIG - 2 - Generic structure of the classical trajectories form- 

metal (light shaded region), ing the basis of the BS scheme. 



This mechanism of 'action cancellation' is the first 
of two basic physical principles determining the phe- 
nomenology of the proximity effect. It not only explains 
the formation of a pair amplitude but also the well known 
sensitivity of the DoS suppression to the classical integra- 
bility properties of the system. This latter point is most 
easily understood within the Bohr-Sommcrfeld approach 
to the problem. The BS approach is based on the con- 
dition that the classical action of any closed trajectory 
traced out by the constituent particles of the system be 
an integer multiple of h. The generic topology of closed 
trajectories of an SN-system is displayed in Fig. 0: for 
sufficiently small energies, e, an electron/hole trajectory 
encompassing two Andreev events closes into itself. The 
total action of such a trajectory is given by 



S 7 (e) - S 7 (-e) 



where the term — h is due to the Andreev scattering phase 
shift iexp(i<f)) [28). Expanding S 7 (e) = 5 7 (0) + T 7 e, 
where T 7 is the time of traversal of 7, one obtains the 
condition 



2T 7 e 



1 

2' 



where n is an integer. This condition implies that for 
small e the structure of the DoS is determined by the dis- 
tribution P(Lry) of classical paths of large length L 1 — 



(See Refs. 12 



J17| , |18| for the precise formula- 
tion of the BS evaluated DoS.) For classically integrable 
and classically chaotic systems, respectively, these distri- 
butions are qualitatively different. More specifically, for 
open chaotic systems, -P(£ 7 ) oc exp(—L 7 /(vFt{)) is expo- 
nentially small (where if is the escape time, see section 
III) whereas for open integrable systems P(L 7 ) decays 
only algebraically. (Here the attribute 'open' means that 
we are dealing with systems that are connected to out- 
side leads.) The exponential rareness of long trajectories 
provides the BS explanation for the smallness of the 
low energy DoS in chaotic systems. 



Notice, however, that irrespective of the structure of 
the system, there will always be some finite probabil- 
ity P(L) of finding trajectories of length ~ fiVp/e. In 
other words, although the BS prediction for the DoS 
may be small, it will never be strictly vanishing (for any 
finite e.) For chaotic systems, this is in contradiction 
with the results obtained by other methods. Specifically, 
for diffusive metallic systems fl|- jio| , pi^j30| , weakly disor- 
dered ballistic systems EaL and 'zero-dimensional' quan- 
tum dot type systems ||12| , [l3f a true DoS gap (strictly 
vanishing DoS for small energies) has been obtained by 
quasi-classical, numerical, and random-matrix methods, 
respectively. 

What is the reason for the discrepancy between these 
results and the BS approach? The answer is that the 
latter does not account for the second key mechanism 
relevant for the physics of the proximity effect, quantum 
diffraction: the discussion above was based on a repre- 
sentation of transition amplitudes in terms of idealized 
classical trajectories. However, in a quantum system the 
uncertainty principle imposes an upper cutoff h on the 
phase space resolution. In practical terms this means 
that the propagation along the trajectories entering the 
BS counting scheme must be interpreted as a caricature 
of the propagation of quantum minimal wave packets. 
The transverse spreading of these packets is at least of 
the order of the Fermi wave length X F (cf. Fig. |). 

Whereas in integrable systems this quantum spreading 
is of no real concern, its consequences in chaotic systems 
are profound. The existence of a globally positive Lia- 
punov exponent A entails that a microscopic initial un- 
certainty of 0(\f ) increases exponentially ~ Aj? exp(iA) 
and, after a time t& = A -1 ln(a/Ai?), reaches the scale a 
of the macroscopic geometry of the system. On time 
scales beyond this so-called Ehrenfest time, a classifi- 
cation scheme based on classical trajectories becomes 
meaningless. The consequences of this spreading of mini- 
mal wave packages have been explored in AL. In that pa- 
per a mechanism was discussed by which the spreading 
leads to the formation of weak localization type quan- 
tum corrections to the classical (Drude) conductance of 
chaotic systems: in the terminology of semiclassics, two 
initially close trajectories split up but eventually recom- 
bine to give an overall phase- insensitive weak localization 
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correction. 

A similar mechanism affects the mean DoS of SN sys- 
tems. However, unlike the weak localization corrections 
to the normal conductance - generally a small correction 
to a large classical background - quantum diffraction in 
chaotic SN systems is of vital relevance for even the lead- 
ing order, qualitative structure of the DoS. In particular, 
it is responsible for the strict vanishing of v{e) below a 
certain threshold energy. 




FIG. 3. Quantum diffraction in a chaotic system: after a 
time set by the largest Liapunov exponent of the system, two 
initially close trajectories diverge. 



For a particular system class, such as diffuse disor- 
dered systems, the effect of these quantum contributions 
can be investigated in quantitative detail. The idea is to 
compare the trajectory oriented picture outlined above 
with quantitative solutions obtained by other methods. 
Indeed, the mean DoS of of a large variety of SN systems 
with diffusive N compounds has been obtained by solving 
the quasiclassical Usadel equations |2^,|6j- |l0| , |32| . Mathe- 
matically, the Usadel equations are second order nonlin- 
ear differential equations for the average Green function 
of the system. They are obtained by a quasiclassical ex- 
pansion of the microscopic Gorkov equations, with self 
consistent account for the (quantum) scattering off im- 
purities. Solving these equations it has been found that 
for systems with an N region of finite extension L, the 
DoS profile universally contains a hard gap of the type 
schematically drawn in Fig. ||. We next address the ques- 
tion of how this phenomenon can be reconciled with the 
semiclassical picture discussed above. 




8* 8 

FIG. 4. Qualitative structure of the mean DoS of a com- 
pact diffusive system in contact with a superconductor. 



For concreteness, let us consider a diffusive system in 
the quantum dot regime, i.e. a regime where the Thouless 
energy E c = hD/L 2 (where D is the diffusion constant) is 
the largest energy scale in the problem. (The discussion 
can be adapted straightforwardly to other types of system 
configurations.) In this case, the Usadel equation is con- 
trolled by the dimensionless parameter n = 7/e, where 
7 is an energy scale representing the coupling strength 
to the superconductor. Although the parameter k need 
not be small, it is nevertheless instructive to expand the 
Usadel equation in k; not to obtain quantitative results, 
but rather to understand the basic physical mechanism 
responsible for the buildup of the gapped DoS. Referring 
to Ref. [|]| for a more detailed discussion we here merely 
sketch the qualitative picture. 

The linearity of k in both the coupling to the super- 
conductor and in 1/e means that each power of this pa- 
rameter corresponds to a single Andreev scattering event, 
coupled to a diffusion mode |33| ] . The diffusion mode de- 
scribes the coherent propagation of a particle and a hole 
through a diffusive background (Fig. ||, right inset), i.e. 
it is the diffusive implementation of the interfering par- 
ticle/hole amplitudes discussed above. A cartoon of the 
correction to the DoS of the lowest order in k is depicted 
in Fig. ||, right. Technically, this correction is described 
by a single diffusive ladder diagram with two Andreev 
scattering vertices [ [l0[ . However, as our present discus- 
sion serves a purely qualitative purpose, we are not go- 
ing to elaborate on the diagrammatic interpretation any 
further. The only thing to be kept in mind is that a per- 
turbative semiclassical/diagrammatic approach to calcu- 
lating the DoS amounts to an expansion in powers of the 
parameter k. 

Already at the next to leading order in k contribu- 
tions of a different topology appear (Fig. [|, center) . The 
distinguishing feature of these corrections is that they 
contain 'junction' points at which a pair of interfering 
particle/hole amplitudes splits up (but later recombincs 
to form an overall phase coherent correction) . These cor- 
rections 

• are not included in a semiclassical approach based 
on a classification in terms of purely classical tra- 
jectories, and 

• are structurally similar to the weak localization cor- 
rections of normal metals. In fact, the basic junc- 
tion region (left inset in Fig. |^) is a real space repre- 
sentation of what in the context of disordered sys- 
tems is called a Hikami box. Within diagrammatic 
approaches these objects appear as vertices in the 
diagrams describing metallic weak localization. 

• Finally, the diffractive splitting of a coherent parti- 
cle/hole amplitude is a process that does not hinge 
on the diffusivity of the medium, but merely on the 
chaoticity of the underlying dynamics. This point 
will be discussed in great detail in section VI below. 
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The main difference to standard weak localization is 
that the DoS 'corrections' depicted in Fig. [| are not 
small. In fact each additional junction introduces two 
more powers of ft, and this parameter can assume arbi- 
trary values. As e is lowered and K /* 1, a situation arises 
where diffraction processes of infinite order have to be 
taken into account, a task which can clearly cannot be 
accomplished within a semiclassical trajectory oriented 
approach. In fact it turns out that for n > 1 expansions 
in terms of trajectories (diagrams) do not even converge 
in principle. 

To understand the origin of the problem, it is instruc- 
tive to compare the semiclassical expansion with the ex- 
plicit solution of the Usadel equation. For the quantum 
dot type system presently under consideration, the qua- 
siclassical expression of the DoS Jl^-[l6|] is given by the 
simple expression, 

v(k) = v Re 7== , 

V 1 — K 

where v$ is the bare metallic density of states (in the ab- 
sence of a superconductor). For our purposes, the most 
important feature of this expression is its BCS type essen- 
tial singularity at k = 1. The structure of this singularity 
implies that (a) n — 1 marks the DoS gap, for (b) values 
k > 1 the DoS strictly vanishes and (c) the gap region 
k > 1 is out of reach of a semiclassical/diagrammatic per- 
turbative expansion centered around the metallic limit 
k — ► 0. For other types of diffusive systems, the rele- 
vant control parameter may be a different one. Also the 
analytical structure of the quasiclassical DoS will in gen- 
eral be way more complicated than the expression above 
[p|-|l0[. However the three points (a)-(c) remain generally 
valid. 

In the rest of the paper we will try to explore to 
what extent these structures generalize from diffusive to 
more general system classes. This analysis will be formu- 
lated within a hybrid formalism comprising quasiclassical 
Green function methods and elements of the supersym- 
metry approach to disordered systems. The relevant the- 
oretical background material is reviewed in the next two 
sections. 








FIG. 5. On the semiclassical interpretation of the solution 
of the Usadel equations. 



III. QUASICLASSICAL FORMULATION 

To promote the discussion of the previous section to 
a more quantitative level, we start out from the Eilen- 
berger equation pi] . While the Eilenberger equation 
may be written in a general form that describes either 
point-like or long-ranged scatterers, we focus in this sec- 
tion on its application to systems with hard boundaries 
and weakly disordered by point-like scatterers ('disor- 
dered billiards'). The case of clean systems with soft 
potentials will be discussed more fully in section VI. 



The general form of the Eilenberger equation may be 
written as 



(ft£-£[a 3 (e + A(r)), ]J g(r, n) = C oii[<7], 

p(r,n)p(r,n) = 1, (1) 

where g(r, n) is the quasiclassical retarded Green func- 
tion of the system |34|]35| ], C the classical Liouville op- 
erator, C = {H, } acting in the phase space of points 
(r,p), and H the Hamiltonian function of the system. 
The coordinates x = (r, n) parameterize &2d—\ dimen- 
sional shell of constant energy Ep = p F /(2m), where 
n = p/pf, Pf is the Fermi momentum and d the dimen- 
sionality of the system. Further, er^ are Pauli matrices 
operating in the two component particle-hole (PH)-space, 
and A = A(cos</kti + sin ^2) is the Nambu representa- 
tion of the phase dependent order parameter. For an 
SN system the support of A is limited to the S region 
where it must, in principle, be determined self consis- 
tently. With regard to the the structure of the proximity 
induced gap in the N region, however, the details of the 
order parameter profile in S are of secondary importance, 
and it suffices to model A by a step function 



A(r) 



A r e S, 
r S N 



For simplicity we will assume that the superconductor 
has a constant order parameter phase = 0. 

The operator O co \\ [g] appearing on the right-hand side 
of equation (|l|) describes the effect of impurity scattering. 
Its functional form depends on the microscopic realiza- 
tion of the scattering potential: 



O co n[g} 



h r dri 

2 J ~Sd~ 
w(n-n') = 2-Kvm\V(n- n') | 2 



[w(n- n').g(r,n'),ff(r,n)], 



where Sd is the surface of the d-dimensional unit sphere. 
Here V is the Fourier transform of the scattering ampli- 
tude and rii is the impurity concentration, while as usual 
the limits of V — > and rn — > 00 are implied, such that 



G 



ni|T^| 2 is constant. The scattering lifetime r is defined 
through w(n — n') by 

1 f dn ' ( <\ 

For the weakly disordered systems considered in this pa- 
per, the scattering time exceeds the typical time of flight 
through the system. 

For finite-size impurities characterized by some Gaus- 
sian potential envelope function, we have 



w(n — n') oc exp 



2 2 

V F a 



n — n 



where a is the typical extent of the scattering region. 
In the remainder of this section we will focus on an 
ultra-quantum limit obtained by sending Ti to infinity at 
fixed Fermi energy Ep and fixed a. (Notice the differ- 
ences with the limit of h — > at fixed electron density 
oc (pF/h) d , usually taken in Boltzmann transport theo- 
ries.) In this limit the form factor w becomes indepen- 
dent of the momentum transfer oc n — n': hence this limit 
corresponds to purely s-wave scattering. The motivation 
for considering the s-wave limit first is that the impurity 
scattering becomes maximally diffractive, making it an 
ideal test case for studying the impact of quantum ef- 
fects on the DoS-profile. Indeed, for s-wave impurities, 
it is the ordinary scattering time r that sets the time 
scale after which classical trajectories are destabilized 
through diffraction and one expects H/t to be relevant 
energy scale for the DoS gap. More complex mechanisms 
of diffraction will be considered in section VI when we 



discuss the proximity effect in chaotic systems. 

As a preparation to the discussion of the full Eilen- 
berger equation (|l|) it is instructive to consider the 
diffractionless limit, for which r — * oo. Physically, this 
limit describes the regime of large energies, e ^> %/t 
(time scales much shorter than the scattering time) for 
which the collision operator is negligibly weak as com- 
pared to both the order parameter and the energy ar- 
gument. In the limit of no diffraction, the Eilenberger 
equation reduces to an ordinary first order differential 
equation that can be solved analytically ]l7j ). Indeed, 
each phase space point x can unambiguously be assigned 
to a classical trajectory x(t) passing through it. The ini- 
tial and final point of this trajectory are set by two points 
of Andreev reflection (cf. Fig. ||). Denoting the time of 
the first (second) Andreev scattering event by — (+)<o/2 
and using that x = {H, x}, the Eilenberger equation in 
the normal region reduces to 

(hdt - ie[a 3 , })g(x(t)) =0, t G [-t Q /2, to/2]. 

In the limit A/e — > oo the solution of this equation is 
given by Q , 

g(x(t)) = sin 9 cos (^J~*J CT 2 + cos #(73 — sin 9 sin (^y^J °' 2 ' 

(2) 



where cos 9 = — itan(eto/7i). This expression contains 
two pieces of information that will be of key importance 
for all what follows: (i) The non-diffractive Eilenberger 
formalism is equivalent to the BS-approach |l2||l3|JT7|Jlq| . 
To derive this equivalence one starts out from the general 
formula for the DoS, 



Ke) = 



Re J dx.tr (g(x, e + )o3), 



where vq is the metallic reference DoS and the phase 
space integral J dx = 1 is normalized to unity. Substi- 
tuting the clean solution for e + = e + iO it is a straight- 
forward matter to show that 



v{e) = 



■KVq 



dx 



00 

E < 

n— — oo 



et (x) 



which is the BS result |IJ,|fJ,|f7|.y . (ii) For energies 
e <C h/to(x) much smaller than the inverse flight time 
through the normal segment of the trajectory, g(x) ps 02 
is approximately 'locked' to the order parameter of the 
adjacent superconductor. This fact will become of impor- 
tance below when we analyze the full Eilenberger equa- 
tion. 

We now proceed to the discussion of diffraction effects. 
For finite r, the Eilenberger equations along different 
trajectories become coupled and no rigorous solution is 
possible. In chaotic systems, however, one may invoke 
statistical arguments to turn the complexity of the prob- 
lem to an advantage. (Notice that under the conditions 
specified above, the presence of impurities alone suffices 
to render the dynamics chaotic on time scales t > r.) 
More specifically, we are going to exploit the features of 
the path length distribution function P(vpto) of systems 
with chaotic dynamics to propose a trial solution of the 
Eilenberger equation that should be close to the unknown 
exact solution. Contrary to the BS result, the DoS de- 
rived from the approximate solution will exhibit a hard 
gap. In the next section we will then show that this gap 
is robust to all orders of perturbation theory around the 
trial solution. 

We begin by noticing that in the diffractive problem 
two intrinsic time scales exist: The scattering time r and 
the 'typical escape time' tf, i.e. the time scale tf = L{/vf 
corresponding to the maximum L{ of the path length dis- 
tribution function P(L). For a ballistic system, we have 
r 3> tf , and so we conclude that three qualitatively differ- 
ent temporal regimes exist: for short times t -C tf, both 
diffraction effects and the differences between integrable 
and chaotic dynamics are negligible, and the DoS is well 
approximated by the BS formula. For intermediate times 
t{ <C t <C r diffraction is still inessential but the struc- 
ture of the function P(L) begin to matter. Specifically, 
the exponential rarity of trajectories with to > tf im- 
plies a strongly suppressed DoS. Finally, for times f>r, 
diffraction is the dominant physical process. 
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Technically, the relevancy of diffraction follows from 
the fact that for small energies e <C Tl/t the scattering 
operator is the dominant term in the Eilenberger equa- 
tion of the normal region. For any given real space point, 
r, the scattering term involves an average over all trajec- 
tories passing through r. Apart from some exponentially 
rare exceptions, the vast majority of these trajectories 
will have a flight time of the order t « tf -C h/e. For 
all these trajectories, the solution of the clean Eilen- 
berger equation is approximately given by 02, i-e. set 
by the bulk value of the order parameter. This justifies 
the approximation of the collision operator by Ocoii ~ 
— ?i(2t) _1 [(T2, ]. Physically, Ocoii then acts like an in- 
duced order parameter of strength « h/(2r). This effec- 
tive order parameter field in the N region describes the 
large angle scattering of Andreev states off the induced 
proximity pairing field amplitude. (Notice the difference 
between a pairing field amplitude, (ipip), and an order pa- 
rameter, which is given by the coupling constant x (ipip)). 

The effect of Ocoii on the solution of the equation de- 
pends on the length of the trajectory under considera- 
tion. For trajectories with t <C r, the clean solution 
(||) is close to the bulk form set by the presence of the 
superconductor at the boundaries anyway, and the extra 
presence of the diffraction operator is of no qualitative 
significance. For the sparse set of exceptionally long tra- 
jectories with to > 2t, however, the boundary condition 
is largely inessential and the equation is governed by a 
competition between the energy and the collision opera- 
tor in the bulk N region. Inspection of the equation shows 
that the solution will approximately be of standard BCS 
type 



?(x) 



e<fi/(2-r) 1 
«o(x)>r ~ ^J{Ti/2t) 2 



iea 3 



•It 



""2 



(3) 



Eq. (9) leads to vanishing DoS for energies e < ?i/(2t). 

Summarizing we have formulated a primitive trial so- 
lution go( x , e) to the Eilenberger equation that predicts 
a gapped DoS. For trajectories io( x ) *C r (energies 
e 3> h/r), go coincides with the diffractionless solution 
(|J). Only for very long trajectories the non-diffractive 
asymptotics is cut off by the collision operator, leading 
to the trial solution of Eq. (||). We emphasize that the 
simple BCS solution (|3|) is specific to the case of purely s- 
wave scattering and becomes unreliable in the vicinity of 
the gap, e « Ti/{2t). In the next section we will explore 
the sensitivity of the result to perturbative corrections 
around the trial solution. It will turn out that the van- 
ishing of the DoS below the gap is a phenomenon robust 
to infinite order in perturbation theory. 



IV. FIELD THEORETIC FORMULATION 

The previous discussion raises the question of whether 
the DoS gap represents a generic feature of disordered 



billiards systems or whether it is an artefact of the BCS 
type low energy trial solution derived above. In princi- 
ple, one may address this issue by using the trial solu- 
tion as a starting point for an infinite order perturbative 
expansion of the Eilenberger equation. (Given that we 
are looking for residual corrections to a vanishing DoS, 
any conclusive perturbative approach has to be driven to 
infinite order.) In practice, however, such a type of per- 
turbative expansion is neither easy to formulate nor does 
it provide much physical insight. Given this situation we 
find it more convenient to follow a different approach. 
This alternative route, which will be based on a varia- 
tional principle, has three advantages: (i) the variational 
approach will provide us with a comparatively transpar- 
ent verification of the robustness of the gap, (ii) it draws 
upon analogies between the curre nt p roblem and more 
general recent developments f^^Jjql regarding the dy- 
namics of chaotic mesoscopic systems, (iii) it will provide 
the theoretical background for the subsequent analysis of 
clean chaotic SN systems. 

What is the meaning of the phrase 'variational prin- 
ciple' in the present context? Within the framework of 
normal systems, Khmelnitskii and Muzykhantskii have 
found [ p5[ a general scheme whereby quasiclassical kinetic 
equations are obtained through a variational analysis of 
certain functional integrals of nonlinear cr-model type. 
Below we will introduce the (straightforward) extension 
of this scheme to the case of superconductor problems. 
Thereafter the functional integral formalism will be ap- 
plied to an analysis of the DoS. By integrating out fluctu- 
ations around the trial solution discussed in the previous 
section, we will find that the gap survives. 

To discuss the variational principle, we need to intro- 
duce some field theory oriented formalism which, at first 
sight, might appear to be quite unrelated to the previ- 
ously discussed material: Consider the functional integral 



Z[J] = J VTe~ s V\ 



where the action S[T] = S Icg [T] + S co n[T] with 

So PI = -tt/5 J dxstrihTasCT- 1 + 

+io- 3 (e + A)Q + J®EflQ\ 
S coU [T} = ^ J ^^ W (n,n')str (Q(r, n)Q(r, n')) . 

(4) 

Here T(x) and Q(x) = T(x.)asT{x.)~ 1 are four dimen- 
sional supermatrix fields in phase space. As before, the 
(unit-normalized) integration extends over a phase space 
shell of constant energy Ep, and the inverse level spac- 
ing at Ep is denoted by p. The rest of the notation 
is adopted to conventions standard in the literature of 
the nonlinear cr-model [B7|: The field T takes values in 
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the (maximally Riemannian subset) of the graded coset 
space GL(2|2)/(GL(1|1) <g> GL(1|1)). In an explicit in- 
dex representation, T(x) = {T aQ ^(x)} where the in- 
dices a, (3 = 1,2 discriminate between the bosonic and 
fcrmiomc blocks of T and the index a, b = 1, 2 is the PH 
index, 'str' is the standard supertrace defined through 
six (A) = — J2 a a A aa .aa(—) a (which differs by an overall 
minus sign from Efetov's convention.) Finally, the two 
dimensional matrix field J — {J a &(x)} serves as a source 
field to generate correlation functions, and Eff is defined 
as (Efl)ij = SuSji. 

An action of the structure of (Q) has first been intro- 
duced in the context of normal systems with weak im- 
purity scattering pq ] . (For an alternative derivation, see 
Ref. ]26[].) For a normal system, A = and the action 
describes the correlated propagation of an retarded and 
an advanced Green function in a ballistically disordered 
environment. (The role of the indices a, b then is to dis- 
criminate between the advanced and the retarded sector 
of the theory.) Technically, the functional integral Z[J] 
is of nonlinear cr-model type, so that Eq. (|]) is commonly 
designated as the action of the 'ballistic cr-model'. 

In the context of superconductivity, the model de- 
scribes a single impurity averaged Gorkov Green func- 
tion. More specifically, 



G(r,n) 



1 



6Z[J] 



.7=0 



tt^O <5J(r, n) 

= (Qn(r,n)) Q , (5) 

where v = p/L d and G(r,n) is the energy-integrated 
Wigner transform of the microscopic Gorkov Green func- 
tion | p5[ . Further, the indices in Qn refer to bf space, 
and (. . .)q denotes functional averaging with respect to 
the source free action S[Q]. From (||) the DoS is obtained 
as 

f=y y*rfrrfnRc(tr (Q u (r, n)cr 3 )) Q . 

The derivation of both the effective action (||) and the key 
relation (|^) is discussed in Appendix [a|. We here merely 
notice that the construction of the model involves essen- 
tial elements of the quasiclassical approach. In fact, the 
derivation essentially amounts to an implementation of 
the quasiclassical construction steps - disorder average, 
Wigner transform, and integration over the modulus of 
the momentum p - into a functional integral approach. 
(The additional graded structure of the functional inte- 
gral is physically of marginal relevance; supersymmetry 
merely ensures unit normalization Z[J = 0] of the source 
free partition function. In fact, this aspect will not be of 
much importance in our further analysis.) 

What is the connection between the field integral ap- 
proach and the Eilenberger equation and, closely related, 
why contemplate a functional integral in a situation that 
is well described by quasiclassics? To address these is- 
sues, we first notice that the functional integral Z[J] can- 
not be evaluated in closed form. It is, however, a good 



candidate for a stationary phase approach. (This is be- 
cause in the present application the coupling constants 
appearing in the action are much larger than unity.) As 
usual with nonlinear cr-models, the action S[T] possesses 
not just one but rather an entire manifold of degenerate 
saddle points (a fact that readily follows from the obser- 
vation that the action is invariant under transformations 
in boson/fermion (bf) space.) We start out by seeking 
for those saddle points Q whose matrix structure is as 
simple as possible. Given the BF invariant structure of 
the action, one anticipates that these saddle points will 
be proportional to unity in bf space, T = {T a bS a p}. 

Defining Q = Ta^T^ 1 , it is a straightforward matter 
to show that the variation of the source free action leads 
to 



5S[T] 

WW 
=> 

where 

Ccoll. 



= 



hLQ(x) - i[a 3 (e + A(x)), Q(x)] = O co u.[Q], 



do! 



to(n,n / )[Q(r, n'),Q(r, n) 



The above equations state that the (bf diagonal) mean 
field configuration of the functional integral arc deter- 
mined by the Eilenberger equation or, equivalently, that 
Q = g is to be identified with the quasiclassical Green 
function. This connection can be made explicit by eval- 
uating (||) on the diagonal saddle point level which gives 

G(r, n) ~ e- s ^Q{r 7 n) = Q(r, n) = g(r, n), 

where the second equality is based on the fact that (due 
to supersymmetry) the action S[Q] = 0. 

Another way of interpreting this result is to say that 
the ballistic cr-model implements a variational principle 
of quasiclassics, a fact first realized in Ref. p^j (with 
reference to normal systems.) In the context of super- 
conductivity this observation is of significant practical 
relevance, the reason being that the variational action 
(|I|) turns out to be a tool more flexibly applicable than 
the Eilenberger equation itself. More specifically, 

• the field integral approach enables one to include 
quantum fluctuations around the mean field ex- 
tremum. As was shown in Ref. |l(J (for the case 
of the dirty limit, A < Ti/t) such fluctuations may 
lead to a renormalization of the leading order qua- 
siclassical results. 

• it is straightforward to extend the field integral for- 
mulation to the computation of the disorder aver- 
age of more than one Green function. This enables 
one to analyze mesoscopic fluctuation effects (!(]] . 
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• in cases where the exact solution of the Eilenberger 
equation is unknown one may still organize the 
functional integral around an approximate mean 
field solution. Provided that this configuration is 
not separated from the true Eilenberger solution 
by points of non-analyticity, no information is lost 
and a (sufficiently accurate) integration over fluctu- 
ations around the trial solution will produce correct 
results. 

For our present application, the third of the above points 
is of relevance. To explore the DoS at low energies e one 
may perform a functional integration around an approx- 
imate solution of the Eilenberger equation. In the next 
section this scheme will be applied to the disordered bil- 
liards discussed in section [II . In section [v| we will then 
extend the discussion to systems with smooth potentials. 



V. SPECTRUM OF AN SN-SYSTEM WITH 
BALLISTIC DISORDER 

In this section we apply the strategy outlined at the 
end of the previous section to explore the DoS of a 
disordered billiard (with s-wave impurities) for energies 
e <C %/t. For these low energies, it is convenient to 
organize the functional integration around the reference 
configuration 



Q r = Rcr 3 R, 1 = cr 2 , 



where the rotation matrix is 



/ 7T 

it = exp yi—ai 



Comparison with the discussion of section |III| (Eqs. (|2|) 
and (||)) indeed shows that for e <C h/r, this configura- 
tion is close to the exact Eilenberger solution (see below 
on what precisely is meant by 'close' in the present con- 
text.) Alternatively, one might organize the integration 
around the trial solution (||). For practical reasons, how- 
ever, Q T = 02 is a more convenient starting point. 

To include fluctuations around Q T , we next introduce 
the parameterization 

Q = RQR-\ Q = Ta 3 T-\ 

where T(x) is a four-dimensional supermatrix field on the 
constant energy shell in phase space. Substitution of this 
ansatz into the action of the functional integral leads to 

SoPI = -TrpJ^y J |±str(?«; F n • TasdrT- 1 - 



Q \rp\ f ^ f dndU ' 



(6) 



str (Q(r,n)Q(r,n') 



The functional expectation value associated with the 
Green function assumes the form 

G(x) =ttv ((RQR- 1 ) 11 (x))q, 



which leads to 

U = ~Y J dxRe(str(<3ii(x)o- 2 ))Q, 

for the DoS in the normal region. In (|J) we have used 
that for the problem under consideration the Liouvillc 
operator C assumes the form -u^n • d T . The symbol J N dr 
indicates that the integration is restricted to the N-region 
of the system (with volume V) . At the boundary to the 
superconductor, we have T = 1, which implies that Q 
becomes tightly locked to its bulk superconductor value, 
(72- Comparison with the general expression (Q) shows 
that (i) due to transformation by R the Pauli matrix 
multiplying the energy argument has rotated to the su- 
perconductor reference point (03 — > —02), and that (ii) 
in Eq. (^j) the appearance of the order parameter is im- 
plicit, namely through the boundary condition imposed 
at the SN interface. 

To prepare the integration over fluctuations, we next 
introduce the standard |57| system of coordinates T — 
1 + W/2, where the matrices W are defined through 



W 



B 



B 



the block decomposition is in ph space and the 2x2 su- 
permatrices B are subject to the (convergence ensuring) 
relation 

B = afBl 

We substitute this representation into the action (0) and 
expand formally according to 

00 

S[T]^S[B,B] = J2S {n) [B,B], 

n=l 

where [B, B] is the contribution of nth order in B, B 
to (||). Similarly, the DoS functional expectation value 
can be represented as 



« 00 

/ dx^(DW[B,B]) Bta , 



(7) 



where the operators [B, B] are defined through the 
series expansion 

00 

str(Qn(x)a2) = ]T£ (2n+1) [£;,5], 



and we have noticed that only contributions of odd power 
in B appear. To prepare the perturbative integration 
over the B degrees of freedom we next consider the second 
order operator [B, B] of the action. It is a straightfor- 
ward matter to verify that the explicit form of [B, B] 
is given by 
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hirp 



dr 
V 



dn dn' 



str 



B(r,n) 



-v F n ■ d r + i J S(n - n') - i ) B(r, n') 



The operator 



—vfu ■ d r + — ) <5(n — n') — — 
t I T 



governing the quadratic action has two key structural 
features (which, as will be discussed in the next section, 
generalize to different system categories): it is (a) real 
and (b) possesses a spectrum that is gapped against zero. 
Whereas the reality of II -1 is manifest, the presence of 
a gap follows from the fact that Dirichlct boundary con- 
ditions B = at the interface to the superconductor 
are imposed. This implies that the phase space zero 
mode -B(x) = const., the only candidate for a zero energy 
eigenmode, does not belong to the representation space 
of n _1 . It is straightforward to verify that all eigen- 
modes with a finite variation have eigenvalues of at least 
0(h/r). 

The features (a) and (b) imply that the operator II 1 
possesses a real and non-singular inverse II (x, x'). As we 
will discuss in a moment, these properties suffice to es- 
sentially determine the structure of the low energy DoS. 

To perturbatively integrate out the -B's, we employ the 
supersymmetry implementation |37 36 of Wick's theo- 
rem: the basic Gaussian pair contraction of two matrices 
B and B is determined by the rules 

(str (B(x)AB(x')A')) = n(x', x) str (A) str (A'), 
(str (B(x)A) str (B(x')A')) = n(x', x) str (AA'), 

where A and A' are arbitrary two-dimensional superma- 
trices, and (...)„ = / V(B, B) exp(-S*( 2 ' [B, B]){. . .) de- 
notes averaging with respect to the Gaussian action. The 
Gaussian average 



(0 ( ( ll) [B,B]0^ 2, [B,B]...) n 

of an arbitrary constellation of rifcth order operators 
can then, at least in principle, be computed as 
the sum of all possible total pair contractions. 

In general, the infinite series emanating from an expan- 
sion of the exponentiated action is far too complex to be 
perturbatively brought under control. With the current 
problem, however, the situation is different. It turns out 
that the series expansion of both the operators entering 
the action and the preexponential terms are sufficiently 
highly structured to make possible some categoric (in- 
finite order in perturbation theory) statements. To be 
more specific, notice that the operators appearing in the 
action fall into two different categories: 

• operators O c [T] = J2 n Ce 2 ™^ [B, B] whose expan- 
sion contains only even powers of B and which 
are real (in the sense that the operator content 



("2 



multiplying the -B's is real). The two represen- 
tatives of this category are the kinetic operator, 
~ str(n ■ Ta^drT" 1 ), and the collision operator 

(<Q)n{Q)x 



str 



• operators of the type i0 o [T] = { 2n+1) [B, B] 
where O is again real but contains only odd pow- 
ers of B. This class of operators is represented by 
the energy term, ~ str((72Q). 

Similarly we notice that the preexponential expansion (Q) 
belongs to the class of purely odd operators. However, 
contrary to the operators appearing in the action, the 
preexponential expansion is purely imaginary (the i's ap- 
pearing in the definition of the Pauli matrix 02.) Finally, 
we repeat that the 'propagator' II is real. 

With all these things in place we now investigate 
the perturbative representation of the DoS expectation 
value. Consider the contribution of any particular oper- 
ator _D( 2n+1 ) [B, B] to the expansion (0). In order to ob- 
tain a non- vanishing contraction, we need to bring down a 
complementary odd operator from the action. The most 
general odd operator that can be constructed from the 
exponentiated action is of the form 

jkQpnx+l) _ _ _ Q(2n k + 1) Q(2 mi ) _ _ _ Q(2m,) ^ 

where I is arbitrary but k must be odd. The total con- 
traction 



D (2n+l)Q(2m+l) 



n (2r lfc + l) n (2m 1 ) 

= i x real 



ol 2mi) 



is then purely imaginary, which implies that no contribu- 
tions to the DoS, Eq. (^), obtain. Summarizing, we have 
found that the mere complex structure of the operators 
coupling to the DoS entails the absence of perturbative 
DoS corrections. Before moving on, let us make a few 
general remarks on the structure of the argument: 

• The construction formulated above explicitly relies 
on the fact that the energy operator in the ro- 
tated representation is purely imaginary. At first 
sight it looks like this makes the derivation explic- 
itly gauge dependent. Indeed, for order parameter 
phases <j> ^ the rotation to the superconductor 
configuration introduces operators into the action 
that are neither purely real nor imaginary. How- 
ever, it is straightforward to verify that (a) the 
phase couples only to the class of odd operators 
and (b) that any non-vanishing Wick contraction 
contains an equal number of phases exp(i^) and 
exp(— i<p). This ensures the gauge invariance of the 
argument. 
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• We emphasize that, in spite of the fact that the con- 
struction is based on the zero energy Eilenberger 
solution g = 02, it is by no means equivalent to 
a mere perturbative expansion of the DoS around 
zero energy. In fact we might have used any other 
low energy approximate solution Q r - in particular, 
the trial solution, Eq. (^), may also have served as 
a reference configuration. As long as Q r is not sep- 
arated from the 'true' Eilenberger solution Q by 
non-analyticities, Q can be reached from Q r by a 
perturbative expansion in W, and the conclusion 
of a vanishing DoS remains valid. The absence of 
non-analyticities for small energies is in turn guar- 
anteed by the presence of the diffraction operator 
O co n. As follows from the structure of the effec- 
tive action, O co \\ causes the induction of a stable 
BCS gap ~ %j (2t) on all classical trajectories in the 
normal region. For subgap energies, e <C h/r, the 
mean field equations of the theory are regular and 
the perturbative scheme works. (The same is true 
for the regime of large energies e^> Ti/t where the 
Bohr- Sommerf eld DoS of the non-diffractive system 
is retrieved.) In the vicinity of the gap region, how- 
ever, the theory is fundamentally unstable and can 
not be used to produce reliable results. 

• The stabilization of the gap by a non-analyticity in 
the functional integral can explicitly be tested for 
the particular case of diffusive SN systems. In the 
diffusive regime, the Eilenberger equation trans- 
mutes to the Usadel equation which, at least for 
a number of relevant geometries, can be solved in 
closed form (cf. the discussion of section ||.) A field 
integral analysis of fluctuations the Usadel mean 
field solutions |l(J then shows that no perturbative 
corrections to the zero subgap DoS exist. The re- 
sponsible mechanism is perfectly analogous to the 
one discussed above. 

We remark also, while we have demonstrated the ab- 
sence of perturbative corrections to the zero subgap DoS, 
this conclusion does not rule out the possibility of con- 
tributions to the subgap DoS of a nonperturbative na- 
ture. For diffusive systems, such contributions to the 
subgap DoS do exist, as pointed out in Ref. (n]] and 
later investigated in detail j38|, and are associated with 
a rare class of states, known as 'tail' states. These states 
are poorly contacted with the superconductor and can 
exist below the gap to give rise to contributions to the 
DoS that are, for diffusive systems, exponentially small 
in the dimcnsionless conductance. In terms of the field 
theory, such states correspond to configurations of the 
Q-matrix that are perturbatively inaccessible from both 
the conventional Eilenberger saddle-point and the metal- 
lic saddle-point of the theory. (More precisely, such states 
can be associated with saddle points that break causality 
in the fermion-fermion sector of the theory.) The study 
of tail states in the ballistic regime, while an interesting 



and unexplored topic, is beyond the scope of the present 
work. 

Finally, to prepare the analysis of soft scattering sys- 
tems, we introduce a bridge between quasiclassics/field 
theory and the qualitative trajectory oriented pictures 
discussed in section ||. In order to keep the discussion 
simple, it is convenient to consider a situation where 
the non-diffractive Eilenberger solution is close to the 
metallic reference point (large energies) and can be in- 
terpreted as the quantum mechanical amplitude of pair 
propagation of ordinary electrons and holes. The exten- 
sion to strongly proximity influenced regimes is techni- 
cally straightforward but less easy to visualize. 

The most basic constituent of the field theory discussed 
above is the propagator II, i.e. the (inverse of the) kernel 
of the quadratic action. Expanding the action around the 
metallic reference configuration Q = 03 it is straightfor- 
ward to show that in the diffractionless limit, the kernel 
is given by 



IT 1 = -j- {~h£ - 2ie) . (8) 



This means that the propagator of the theory in phase 
space, II(x, x') is just the e-Fourier transform of the clas- 
sical time evolution operator, 



9 r°° 

n(x, x') = -4- / dte 2Kt ' n (e ct ) (x, x') = 

■Khp J 

9 r°° 

= / dte 2fc */**(x-x / (t)). 
■Khp J Q 



By expanding the functional integral in terms of the B 
field, one may now describe the effect of Andreev scat- 
tering and/or quantum diffraction in terms of perturba- 
tive superpositions of classical propagators. For exam- 
ple, expanding the functional to lowest order around the 
quadratic kernel the first correction to the metallic DoS 
is found to be 



Sp = -npA 2 i>o Im — / dxirfx 2 II(xi, x 2 ), 
oe Js 



where the II describes the propagation along the classi- 
cal trajectory of an Andreev bound state (Fig. |[ right), 
and J s denotes an integration over the superconductor 
interface. 



12 




FIG. 6. First order quantum correction to the classical 
propagator in a SN-system with hard impurity scattering. In- 
set: solid line: propagation amplitude of an electron, dashed 
line: amplitude of the Andreev reflected hole state, dot: 
s-wave impurity scatterer. 

In the same manner one verifies that, at next order in 
an expansion of the functional integral in the generators 
B, a quantum correction to the classical propagation of 
the type 



<ttl(xi,x 2 ) = -2- 



hA 2 (irp) 3 f dr f dridn 



T Jn V J £>d J s 
xll(r, n; xi)n(x 3 ; r, n)n(r, n'; x 4 )n(x 2 ; r, n') + . . . , (9) 



dx-idx.4 x 



appears. (Here the ellipses denote contributions of sim- 
ilar architecture but somewhat different structure in 
phase space.) 

The integration J s extends over (the phase space 
points on) the superconductor interface. A cartoon of 
this process is depicted in Fig. |[ where the inset repre- 
sents a visualization of amplitudes involved in the diffrac- 
tive s-wave scattering. We notice that SH represents 
the quantitative formulation (for a hard scattering ballis- 
tic environment) of the diffraction corrections discussed 
qualitatively in section [n| Through such processes, an 
exceptionally long classical trajectory (quantitatively de- 
scribed by the propagator n(x, x')) establishes contact 
with the superconductor on time scales r. Substitution 
of the modified propagator into the formula for 5v gives 
the first quantum correction to the DoS. 

Technically, the quantum contribution to II can be 
interpreted as some 'self energy' correction of strength 
~ Ti/t. On small energy scales e < h/r the action of the 
diffraction operator is no longer perturbatively accessi- 
ble. Nonetheless, the concept of a low order expansion 
around the limit of classical propagation remains a use- 
ful concept: In cases where the mechanism of quantum 
diffraction is less trivial than mere s-wave scattering, a 
perturbative analysis can be used to determine the char- 
acteristic energy scale beyond which the quantum propa- 
gation becomes dominated by diffractive mixing between 
classically isolated regions in phase space. The general 
structure of the theory implies that this time scale sets 



the width of the gap. We will come back to this point 
in the next section, when the soft scattering problem is 
discussed. 



VI. CLEAN CHAOTIC SYSTEMS 

In this section we turn our attention to the proximity 
effect in clean chaotic systems. As before we focus on 
billiard type systems in which the average time of flight 
if represents the shortest time scale in the problem. Wc 
assume that the scattering potential is smooth. (In bil- 
liards, i.e. systems where the characteristic length scale 
of the scattering potential a coincides with the system 
size, 'smooth' means a ^ Xf- For a definition applicable 
under more general conditions, see |p4j| .) 

The traditional form of the Eilenberger equation, 
Eq. (|l|), applies to billiard type systems subject to weak 
impurity scattering. Its naive generalization to generic 
clean systems characterized by a set of phase space coor- 
dinates {x} is given by 



h£-i[a 3 (e + A(x)), ] U(x) = 0, 



(10) 



with the standard non-linear normalization <?(x)<?(x) = 
1. This equation differs from ([!]) by (i) the absence of 
an impurity scattering operator O co \\ and (ii) the gener- 
alization from i^n- d r to an arbitrary Liouville operator. 

A moment's thought, however, shows that Eq. ( fTo| ) 
poses a serious problem: as discussed in section III, the 
diffractionless Eilenberger equation reproduces the BS re- 
sult for the DoS. (No specific assumptions on the struc- 
ture of the Liouville operator were made in that argu- 
ment.) Relying on (|To|) one would thus have to conclude 
that the proximity DoS of clean chaotic systems is not 
gapped, i.e. that the DoS gap must be a peculiarity of 
disordered systems. This hypothesis, however, can be 
dismissed on the basis of a simple argument: being a pos- 
itive definite quantity, the DoS of disordered systems can 
only vanish on average, if it vanishes for each individual 
realization. However, a given realization of a disordered 
system is nothing more than a Hamiltonian chaotic sys- 
tem. From this argument one expects that the proximity 
gap is a generic signature of chaotic dynamics (save, per- 
haps, for chaotic systems - unknown to us - that do not 
fall into the universality class of standard Hamiltonian 
systems with chaotic potential scattering.) 

The resolution of the paradox is that for clean chaotic 
systems the derivation of ( [To|) is ill-founded. Indeed, 
the derivation of the Eilenberger equation is based on 
the presumption that the quasiclassical Green functions 
varies slowly on the scale of the Fermi wave length. For 
disordered systems, the ensemble averaging over realiza- 
tions makes sure that this requirement is met. With 
chaotic systems, however, the situation is different. The 
phase space eigenfunctions of the chaotic Liouville oper- 
ator are singular and fluctuate on all length scales. This 
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invalidates any semiclassical gradient expansion, unless a 
smoothing regularization mechanism is active. 

Before discussing the issue of regularization further, we 
switch back to the language of the field integral formula- 
tion. The variational action corresponding to the clean 
Eilenberger equation is given by So [T] of Eq. (||) . After 
the rotation T — > RT, R = exp(i7rcri/4) to the super- 
conductor reference point {Q = 02), the action simplifies 
to 

S [T] = -np J dxstiihT^CT- 1 - iea 2 Q^ (11) 

where coupling to the superconductor is implicit, through 
the boundary condition T = 1 at the SN-interface. 

Within the context of the general ballistic er-model, the 
necessity to regularize the effective action So has been 
pointed out earlier ||^,^,^]. More specifically, it has 
been noticed that the phase space singularity alluded to 
above invalidates the gradient expansion on which the 
derivation of the effective action So is based. To give the 
effective action some meaning, a regularizing contribu- 
tion SVeg has to be added. Frustratingly, however, no gen- 
erally accepted consensus regarding the concrete imple- 
mentation of a regularization operator has been reached 
as yet. In fact, various different schemes have been pro- 
posed in the literature: in the early work pt| it was sug- 
gested to regularize the action by adding a second order 
differential operator of infinitesimal coupling strength. 
In Ref. |Q it was argued that a Hamiltonian chaotic 
system will se//-generate some regularization if only the 
mechanism of quantum dispersion is consistently taken 
into account. In lack of a fully microscopic description 
of diffraction within a quasiclassical formalism, it was 
suggested to model phenomenologically the effect of the 
quantum spreading of minimal wave packets. This was 
done by adding the diffusion type second order operator 

SrM = —£rl dxstr(d n Qd n Q) 

to the effective action, where the value of the character- 
istic scattering time r r = X 2 vf/^f and A is the largest 
Liapunuov exponent of the system (i.e. unlike in Ref. 
pfjj, the coupling strength of the regulator was kept fi- 
nite.) Physically, the operator S rcg [Q] describes the ef- 
fect of some soft (low momentum transfer) scattering off a 
stochastic potential; this may be seen readily by compar- 
ison with the collision term S co \i in the action (Q), in the 
case where w(n—n') is strongly peaked around n— n' ~ 0. 
We re-emphasize that this potential is fictitious and only 
serves to model the spreading inherent to the propagation 
of any quantum minimal wave package. In contrast, Ref. 
p9| argued that a chaotic system will not auto-regularize 
itself and that some external smoothing is needed. It was 
suggested to effect this regularization by adding to the 
action an elliptic second order differential operator cou- 
pled to some weak stochastic background field. Again in 



contradiction to Ref. |24j, it was conjectured that a cou- 
pling strength vanishing in the classical limit would be 
sufficient to generate regularization. 

Although no general consensus regarding the practical 
realization has been reached, it is clear that the ballis- 
tic cr-model needs regularization through some positive 
second order differential operator (or any other suitable 
scheme.) In view of this situation, we here take a prag- 
matic point of view and add to the action So the contri- 
bution 

SVeg = -^fBl J dxstr (diQdjQ), (12) 

where Dij is a positive matrix whose characteristics will 
be left unspecified for the moment. 

Before moving on it is perhaps worthwhile to re- 
emphasize that this regularization problem is in no way 
particular to the field theoretical approach. It rather af- 
fects all semiclassical expansion schemes, including those 
underlying the derivation of the Eilenberger equation. In 
applications involving disorder, the averaging over a con- 
tinuous set of potential variables effectively smoothens 
out all singularities, and the derivation of effective low 
energy theories is not met with difficulties. For individual 
representatives of chaotic systems, however, the problem 
is virulent and regularization must be an integral part of 
any quasiclassical approach. We re-emphasize that the 
actual mechanism by which regularization is effected 'in 
reality' is not well understood at present. 

We now turn to the discussion of the effective action 
S = S + Sreg determined by Eqs. ([H]) and flU), respec- 
tively. Comparison with (gj) shows that S bears much 
structural similarity with the action of the previously 
discussed s-wave scattering system. The essential differ- 
ence lies in the structure of the scattering operators S co ii 
and S rog , respectively. Whilst S co \\ describes large mo- 
mentum transfer scattering at a rate determined by 1/r, 
the diffractive scattering modelled by S rcg is soft and of 
a strength determined by the matrix D. Naively, one 
might conjecture that this difference in the microscopic 
structure cannot be of much concern and that S reg acts 
qualitatively similar to S co n. This anticipation is wrong. 
In contrast to S co n, the operator S reg leads to a forma- 
tion of a gap only for systems with underlying chaotic 
dynamics. Further, the width of that gap depends only 
logarithmically on the coupling strength, D. 

To understand the action of the diffraction operator 
we consider the soft scattering analogue of the perturba- 
tive picture developed in the end of the previous section. 
Expanding the action to second order in the generators 
B, it is straightforward to verify that, in the presence of 
soft diffraction, the classical propagator (||) generalizes 
to 

n^ 1 = ^ (-7LC - 2ie - HDijdidj) . 

One further verifies that the first order quantum cor- 
rection to the propagation along a classical trajectory is 
given by 
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<m(xi,x 2 ) = 8{np) 3 A 2 hD l:j J dx J dx 3 dx 4 x (13) 

x<9 Xi n(x, xi)n(x 3 , x)n(x, x 4 )9 Xj n(x 2 , x), 

(compare with expression (^) for <5IT in the case of 
hard, s-wave scattering). At first sight, the structure 
of this expression suggests that upon sending the regula- 
tor strength to zero, the quantum corrections to classical 
propagation vanish linearly with D. However, this antic- 
ipation is too naive; it does not account for potentially 
singular behavior of the x-derivatives acting on the prop- 
agator n. 

To better understand the behavior of the quantum cor- 
rection for small D, we need to explore the structure of 
the regularized propagator II. Consider the time repre- 
sentation of the defining equation of II, 



np 



(d t - £ x - AAAJn^x'; t) = 5{x - x')5(t), 



(14) 

where t > and the subscript in £ x indicates that 
the Liouvillian acts on the argument x. For notational 
convenience, the factor of two multiplying the energy 
argument e has been included in the definition of the 
time Fourier argument t. In the absence of diffraction, 
n D=a (x,x';t) = (2/tt P )5(x - x'(t)). For finite D tJ , U 
is no longer purely deterministic but rather describes a 
stochastic process of diffusion centered around the clas- 
sical trajectory x'(t). For sufficiently small times, the ar- 
gument x is close to x'(t) and the equation for II can be 
linearized in y = x — x'(t). Introducing C = Iijd Xj Hd Xi , 
where the matrix 



a 



-1 



made: Eq. (|l5|) is governed by a competition between 
the diffusion operator D i jdidj 1 and the drift operator 
\1, J),, //,. By definition, the time averaged eigenvalues 
of the 'drift matrix' M are the Liapunov exponents of 
the chaotic scattering problem. On short time scales, 
the dynamics of the equation is controlled by the diffu- 
sion term and we expect the distribution II to describe 
a stochastic diffusive spreading around the deterministic 
solution x'. After a certain threshold time, a crossover to 
a regime controlled by the drift operator takes place. On 
larger time scales, the distribution II describes some ex- 
ponential separation away from the trajectory x', where 
the rate of separation is determined by the Liapunov ex- 
ponents. 

To formulate these concepts in a somewhat more quan- 
titative manner we make the assumption (acceptable on 
short time scales) that the matrix M is approximately 
static. (In fact, all what we really need to require is 
that the largest Liapunov exponent, A, and the associ- 
ated direction in phase space are constant.) For a time 
independent matrix M, Eq. ( |l5| ) can be solved in closed 
form Q: 

n (y.*) = i0(27r)( t + i)/ 2 ( de ^ft))' 1/2e ' |yr<T(t) " ly ' 

where the matrix-kernel a is given by 

a (t) = 2 ]T —r— - V (v" T DVV T . 



at) 



A Q + A/3 



Here v Q and A Q are eigenvectors and eigenvalues of 
the drift matrix M, respectively. The enumeration of 
the eigenvalues is ordered such that A = Ao > Ai > . . . > 
A2d-i- Finally, re-installing our original set of variables, 
we arrive at 



implements the symplectic structure of the phase space, 
the linearized equation for II assumes the form 

(d t - £ x '( t) - M,,i),, a, - D. J),, (),, ) II(x'(i) + y,x';t) 

= -S(y)S(t), 

np 



with Mi, 



HJ = Iikdl ]Xk H\^, {t) . Noticing that for any 

function f(x'(t)) that is not explicitly dependent on t, 
we have (dt — £ x /(t))/(x'(t)) = 0, the differential equa- 
tion can be rewritten as 



(d t ~ Mijyjdy 



n 8 8 



n(y;i) 



7T/9 



-S(y)6(t), (15) 



where the dependence on both x'(t) and the initial con- 
dition are implicit, and the time derivative acts only on 
the explicit time argument in II(y;t). Eq. ( |l5|) is a dif- 
ferential equation of Fokker-Planck type. Owing to the 
complex time dependence of the matrix M it is, of course, 
not analytically solvable. Nonetheless, some qualitative 
statements regarding the structure of the solution can be 



n(x,x',i) = p(27:)(2d+1)/2 (deta(t))-^ x 



xe -i(x-x'(t)) T ^(t)" 1 (x-x'(t)) 



for the short time approximation to the propagator LT. 
From this expression we recover the temporal picture dis- 
cussed qualitatively above: For times shorter than X , 
d ~ Dt and the solution describes diffusion around the 
classical trajectory x'. On time scales larger than A -1 , 
some of the matrix elements of a begin to diverge expo- 
nentially. Defining xq = v 0T x as the projection of x onto 
the direction of strongest exponential divergence, we find 



n 

A 



V 4 , 



(16) 



where D a = v Dv u denotes the relevant element of 
the diffusion matrix. Thus, for t > A -1 the propagation 
is no longer confined to the vicinity of x'(t) but rather 
branches out exponentially, at a rate determined by the 
largest Liapunov exponent of the system. We emphasize 
that this type of behavior is specific to chaotic systems. 



15 



For an integrable system, no globally positive exponents 
driving the separation process exist, and, the weak dif- 
fusive smearing left aside, II does not deviate from the 
predestined trajectory x'. 

We now turn back to our discussion of the quantum 
correction, STl of Eq. (|l3|). Specifically we wish to un- 
derstand how much time it takes before an exceptionally 
long classical trajectory diffractively spreads out into the 
phase space regions strongly coupled to the superconduc- 
tor. The first thing to be observed is that the prefacto- 
rial dependence of 511 in Eq. (|l^) on the diffusion matrix 
Dij is largely cancelled by the derivatives acting on the 
propagators II. Indeed, the diffraction correction at fixed 
energy, e, is determined by the time Fourier transform of 
the propagators II(x, x'; t). Parametrically, the contribu- 
tion of short times t < A^ 1 to the time Fourier integral 
of <9 X IT is given by 

9 x n ~ (x - xO^n ~ ff-^n ^ (Dty^u, 

implying that the /^-dependence introduced by the two 
derivatives cancels against that of the prefactor. From 
this rough estimate we conclude that it will not be the 
naive timescale related to the coupling strength of the 
operator that sets the effective rate of diffraction, but 
rather some other time scale. 

To determine that scale we apply a modified version of 
an argument first formulated by AL within the context 
of normal metal weak localization: let us focus attention 
on the coordinate xq of strongest exponential divergence. 
For the sake of concreteness, we make the (inessential) 
assumption that Xq represents some real space coordi- 
nate. From the solution for II we find that for short time 
scales t < A -1 , ((xq — x' (t)) 2 ) w Dot. On the other 
hand, the dispersive quantum spreading of a Gaussian 
minimal wave packet centered around the classical tra- 
jectory would be estimated by ((xq — x' (t)) 2 ) ~ vpXpt. 
We now follow the general philosophy of AL and postu- 
late that the mere function of the regularization operator 
is to phenomenologically mimic this effect. This leads to 
the estimate D Q ~ vpXp pl[ . 

Having parametrically estimated the strength of the 
diffraction operator, we next consider the quantum cor- 
rection, SU of Eq. (p"3|), in the 'Liapunov regime' of time 
scales t > A -1 . To couple a formerly isolated excep- 
tionally long trajectory x' : Xi — > X2 to the background 
of average trajectories and/or the superconductor inter- 
face, the propagators II(x3,x) and II(x, X4) 'branching 
out' from the point x of diffractive scattering, have to 
leave the immediate phase space vicinity of the refer- 
ence trajectory. A conservative estimate of the threshold- 
deviation necessary to initiate that coupling reads as 

<(*o-^)) 2 )«fVW, 

i.e. we demand that the spreading reaches a scale compa- 
rable with the system geometry. Finally, estimating the 



Liapunov exponent as A ~ vp/a, and using the above 
expression for D 0l we arrive at 

-Hi)* 

which is the celebrated Ehrenfest time fg. 

Summarizing we have found that for chaotic systems 
with smooth potentials, £e is the time scale at which 
quantum corrections begin to significantly affect the clas- 
sical propagation along ballistic trajectories. On time 
scales larger than £e any formerly isolated long classical 
trajectory 'merges' into the background of surrounding 
trajectories of average flight time U (where tf is the es- 
cape time). According to the general structure of the 
theory, the inverse of this threshold time, h/t^, sets the 
width of the DoS gap. 

It is instructive to study this mechanism more explic- 
itly on the basis of the Eilenberger equation. This is 
done in appendix [b|, where we consider the prototypical 
example of an infinitely long trajectory x' that, classi- 
cally, would never hit the superconductor interface. In 
the absence of diffraction, the Green function along that 
trajectory is metallic, g(x',e) — 03, and gives a finite 
contribution to the DoS. For O co \\ weak but finite, the 
behavior of g close to x' is described by a Fokker-Planck 
equation similar to the one discussed above. The so- 
lution of this equation, detailed in Appendix shows 
that for energies e ~ Ti/tE, a crossover from 03 to the 
off-diagonal configuration set by the surrounding super- 
conductor takes place. For still smaller energies the DoS 
vanishes. 

Before concluding, let us make some general critical 
remarks on the status of the results obtained in this sec- 
tion. Unlike in the case of hard point-like scatterers, the 
derivation and solution of quasiclassical equations (alias 
the mean field equation of the ballistic a-model) relied 
on various phenomenological assumptions: 

• A number of assumptions regarding the nature of 
the Liapunov exponents of the system, and the so- 
lution of the linearized equations of motion have 
been made. Furthermore, some relevant classical 
and quantum system parameters have been esti- 
mated in a crude manner. 

• The all-important diffraction operator has not been 
microscopically derived but was introduced on phe- 
nomenological grounds. 

We believe that the first of these points is not serious. 
The reason is that the characteristic energy scale h/t-E 
depends only logarithmically on the modelling of the dy- 
namics, and on the values assigned to the parameters 
D, A, a, etc. As long as a globally positive Liapunov ex- 
ponent exists, an energy scale of the order of h/t-E will 
be generated. The second point, however, is more im- 
portant. The absence of a consistent microscopic deriva- 
tion of a diffraction operator indeed represents a defi- 
ciency of the current approach and, for that matter, of 
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the general field theory approach to quantum chaos. On 
the other hand, the presence of a diffraction operator, of 
the coupling strength specified above, is required for ba- 
sic physical reasons. We therefore believe that, at least 
phenomenologically, the current approach stands on solid 
ground. 

Finally, we notice that, unlike in the case of the hard 
scattering problem, it was not possible to bring the Eilcn- 
berger equation into a form manifestly equivalent to a 
bulk BCS equation with induced order parameter. This 
implies that the current approach cannot rigorously es- 
tablish the existence of a hard gap edge for the soft scat- 
tering problem. What we could see is that on energy 
scales of the order of the inverse of the Ehrenfest time, 
a quantum mechanism suppressing the DoS becomes op- 
erative. The fact that for lower energies the solutions 
of the Eilenberger equation undergo a rotation towards 
the bulk superconductor configuration (cf. Appendix [B|) 
implies that the BS approach qualitatively overestimates 
the low energy DoS of chaotic systems. 



VII. SUMMARY 

In this paper we have discussed various physical mech- 
anisms relevant to the phenomenology of the proxim- 
ity effect in chaotic SN-systems. Considering systems 
where chaotic dynamics is induced by a weak amount 
of point-like impurities as a test case, we have demon- 
strated that (a) the existence of a hard gap in the DoS 
is directly related to an essential singularity in the qua- 
siclassical equations for the Green function and (b) that 
these singularities can be interpreted semiclassically as a 
non-perturbative accumulation of quantum interference 
processes akin to weak localization corrections in normal 
metals. These contributions are effectively summed up 
in the solutions of the quasiclassical equations. Although 
closed solutions can be obtained only in exceptional cases 
(for example in the excessively investigated case of sys- 
tems with diffusive dynamics) we have shown that the 
strict vanishing of the DoS for low energies follows from 
combining elements of quasiclassics with the formalism 
of the ballistic cr-model. While all perturbative correc- 
tions to the DoS vanish below the gap, we have remarked 
that exponentially small contributions due to tail states 
are not ruled out by our analysis, but remain beyond the 
scope of this work. 

Turning to the far more complex case of clean chaotic 
systems we have discussed various difficulties obstructing 
the construction of quantitative analytical approaches. 
The perhaps most fundamental problem is that the as- 
sumption of a spatially continuous Green function under- 
lying, for example, the quasiclassical formalism fails (due 
to the singular distribution of the phase space eigenhmc- 
tions of the Liouville operator.) In principle one may 
evade this difficulty by averaging the theory over some 
set of external parameters. Structure and strength of 



a minimum set of averaging parameters needed to en- 
sure continuity are not known at present. However, tak- 
ing a pragmatic point of view and assuming continuity 
of the Green function, it can be demonstrated that for 
energy scales smaller than the inverse of the Ehrenfest 
time, quantum diffraction mechanisms become operative 
pi| . For those low energies the DoS of SN systems is 
suppressed below the value predicted by the semiclassi- 
cal Bohr-Sommerfeld approach. (Here it is assumed that 
the Ehrenfest time is larger than the average classical es- 
cape time from the normal metal region. The opposite 
regime has not been investigated so far.) Most likely, the 
inverse of the Ehrenfest time marks the position of a hard 
gap edge although, due to the problem stated above and 
other difficulties related to the formulation of quantita- 
tive theories of chaotic systems (see the text), this could 
not be rigorously shown. In consideration of the diffi- 
culties mentioned above, and the fact that established 
semiclassical techniques based on counting classical tra- 
jectories fall short of producing a hard gap, a quantita- 
tive theory of the proximity effect remains a challenge for 
future approaches to chaotic sytems. 
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APPENDIX A: DERIVATION OF THE 
BALLISTIC SIGMA-MODEL 

In this section we derive the effective action (||) of a 
weakly disordered sample in the presence of a supercon- 
ducting order parameter. While, as discussed in the main 
text, an action of this type for a normal metal was first 
proposed by Muzykantskii and Khmelnitskii [po| , we fol- 
low here more closely the derivation of Andreev et al. [ p6[ , 
with suitable modifications: we allow for the presence of 
a superconducting order parameter by introducing parti- 
cle and hole components in the functional integral. The 
method relies on an averaging over a set N 3> 1 of energy 
levels around the Fermi surface to construct the cr-model 
in the absence of disorder. The derivation employs fur- 
ther procedures that are central also to the quasiclassical 
approximation, including a Wigner-transform and inte- 
gration over the momentum coordinate normal to the 
energy shell. The disorder is treated as a perturbation to 
derive a term resembling a collision integral in the action. 
Since we present the generating functional for a single 
retarded Green function, the role of the usual advanced- 
retarded components in the field integral are taken over 
by particle- hole components. For simplicity we restrict 
ourselves to the case of broken time-reversal symmetry, 
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although this simplification is not essential and the gen- 
eralization to other symmetry classes is straightforward. 
The starting point is the Gorkov Hamiltonian, 



H = H + V(r) + Aa 3 ; H = 



2m' 



The generating function for the single-particle Green 
function may be expressed as a field integral in the form 



x 



x exp 



J dY^{v)L{G- 1 {e )+iJEf l )^{i 



where G~ 1 (eo) = e — ecr 3 — H. ^ represents a 
four-component superfield, containing particle-hole and 
boson-fermion sectors. The energy e represents an en- 
ergy difference between particle and hole sectors. The 
matrix L is defined as 



L 



1 



with the block decomposition being in particle-hole 
space. Two dimensional matrix field J = { J a fc(x)} serves 
as a source field to generate correlation functions, and 
Efl is defined as (Efl)ij — SuSji. Next we employ Gaus- 
sian averaging over energy levels according to 



dec 



(2 7 riV 2 ) 1 /2 



exp[-(e -6 F ) 2 /(27V 2 )] 



This induces a nonlocal interaction term in the generat- 
ing functional: 



(Z(J)h 



Sn 



L»*exp 



J dr^(r)Lx 

x(G- 1 (e F )+iJEf 1 )^(r) - 5 int [*] 

N 2 / f t N 2 

dr*t(r)L*(r) 



The interaction term is then decoupled with a matrix 
Q(r,r') that is nonlocal in space: after integrating over 
we find 



(Z(J))e 



DQexp 



^str r Q 2 



str r ln (Gq 1 {Q) - a 3 (e + A) - V + UEfl 



(Al) 



where Gq 1 = tp — Hq — iNQ and Q obeys the symmetry 
= LQL, where the dagger denotes a Hermitian con- 
jugate. The strategy now is to simplify the action so as 
to isolate fluctuations associated with energies much less 
than the Fermi energy. To do so we search for a value 
of Q that minimizes the action and expand in fluctua- 
tions around this saddle-point. In addition, we consider 



first the totally clean case (V(r) = 0), while later we in- 
clude the effects of ballistic disorder. The saddle-point 
equation now reads as 



QqG (Qo) 
with the explicit solution 



. £f — Ho 
' 2N 



1 



-iN, 



cf — H 
2N 



1/2 



03- 



To describe the low-energy, lon g- w avelength fluctuations 
of Q around the saddle-point |42| . we may write Q — 
TQqT" 1 , where T obeys the symmetry T' = LTL, and 
expand the action to leading order in l/N: we find 



S[Q] = --str 



Q (a 3 (e + A) - iJEH - T^H^T] 



(A2) 



At energies far below the Fermi energy, the leading con- 
tribution to the above functional integral may be de- 
scribed by the semiclassical approximation (note that 
the use of the term 'semiclassical' has an entirely dif- 
ferent meaning from that used in the BS approach; here 
it stands for an expansion of operator traces in powers of 
h). To define this approximation, we apply the Wigner 
transformation 



C(R,p) = J d(6r)e ip - 5r / h O 



R 



5r 



,R 



y 



According to the semiclassical approximation, the 
Wigner transform of a product of operators may then 
be approximated as (0i02)(r, p) = Oi(r, p)C 2 (i", p) + 
fr{Gi(r, p), C?2( r ) p)} + . . ., where {•, •} denotes the clas- 
sical Poisson bracket. Specifically, the semiclassical ap- 
proximation of the commutator with the Hamiltonian 
reads as 

[H ,T] -> -ihCT = -ih{T,H }, 

where C is the Liouville operator. Further, assuming 
that (a) fluctuations transverse to the constant energy 
shell Hn — const, are frozen out and (b) Q Q w a 3 (see 
Ref. ^tf^ for a discussion of these points) we may write 
the semiclassical approximation of the matrices Q as 

Q(x) =T(xy 3 T(x)-\ 

where x is a shorthand notation for the 2c?— 1-dimcnsional 
set of phase space coordinates of the constant energy 
shell. For a billiard type system with hard boundary 
scattering x = (r, n). 

Performing the integration over the coordinate trans- 
verse to the energy shell (for constant Q), we obtain from 
Eq. (A2) the So contribution to the action contained in 
Eq. ( This part of the action is applicable to a clean 
sample in the presence of energy averaging. 
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To include the effects of ballistic disorder, we may now 
reintroduce the impurity potential V(r). Although for 
short range impurities this part of the Hamiltonian may 
not be treated by the semiclassical approximation, in the 
limit of weak disorder it may be treated as a small per- 
turbation in a similar manner to the symmetry break- 
ing terms JEfl and eos appearing in the action (Al). 
Expanding the logarithm to second order in V and av- 
eraging over impurity configurations then leads to the 
collision term Scoii in Eq. 

It may be verified (see main text) that the saddle-point 
equation for Q that minimizes the effective action (Q) 
is precisely the the quasiclassical Eilenberger equation, 
Eq. (Q). The strength of the field theory approach how- 
ever is that it enables evaluation of the single-particle, 
impurity-averaged Gorkov Green function even beyond 
the quasiclassical approximation. Differentiation of the 
partition function with respect to the source term, J, 
leads to 



G(r,n) = 



1 



dZ[J] 



two dJ(r,n) 



,7=0 



= / DQe 



-S[Q] 



<2n(r, n), 



(A3) 



where G(r, n) is the Wigner-transformed and energy- 
integrated form of the Gorkov Green function, and the 
11 indices refer to bf space. While evaluation of the in- 
tegral in Eq. ( |A3| ) at the saddle-point level reproduces 
the quasiclassical approximation for G(r, n) ~ Q(r,n), 
integration over the entire Q-manifold in principle repro- 
duces the exact expression for G(r, n). 



APPENDIX B: APPEARANCE OF THE 
EHRENFEST TIME WITHIN THE 
QUASICLASSICAL FORMULATION 

In this appendix we explore the influence of the colli- 
sion operator on the solution of the Eilenberger equation 
in a chaotic environment. Specifically, we consider the 
equation in the vicinity of a classical trajectory x' of in- 
finite length, i.e. a trajectory that never hits the super- 
conductor interface. This trajectory is meant to serve as 
a caricature of the long trajectories that support the BS 
density of states at low energies. 

As long as the collision operator is not included, the 
quasiclassical Green function on our trajectory assumes 
its metallic value, g(x',e) = 03. For coordinates x away 
from the reference trajectory and/or in the presence of a 
finite collision term we expect deviations from 03. To pa- 
rameterize these configurations we introduce the param- 
eterization g(x) = e w ( x 'aae~ w ( x ' , where W is a 2 x 2 
antiunitary matrix generating rotations around 03. For 
all what follows it will be convenient to parameterize W 
as W — W + + W- , where 



W± = iw±-^=(<r 1 =F icr 2 ), 



w± e K. 



Anticipating that in the vicinity of x', W will be small, 
we next expand g as g = 03 + 2Wa 3 + . . . and linearize 
the Eilenberger equation 

(hC - ie[a 3 , ])g(x) - TiD^di (g(x)djg(x)) = 0, 

in W. The result reads as 

(±h£ + 2ie + hi ),,»;,) w± = 0. 

To make progress with this equation, it is again conve- 
nient to transform to time space. Fourier transforming 
we obtain 



(±C - dt 



Dijdfj) w± 



0, 



where, as in the main text, the factor of 2 multiplying the 
energy argument e has been included in the definition of 
t. To give this equation some meaning, initial conditions 
have to be specified. To this end we introduce a coor- 
dinate system where the first of the 2/ phase space co- 
ordinates y, yo, yi, ... , j/2d-i parameterizes our reference 
trajectory, whereas the remaining coordinates span the 
complementary region. We next assume than all possible 
choices of initial conditions w±(x, t = 0) = w±fl(x) have 
the common feature that w±^q increases up to w± = 0(1) 
as any of the coordinates xq, . . . ,X2d-i approaches the 
characteristic scale of the potential, a. What this means 
is that for phase space points farther away from the ex- 
ceptional trajectory than a, the Green function is locked 
to the superconductor order parameter. We notice that 
neither the specific choice of the threshold value xt m a 
nor the rigorous implementation of the initial condition 
sensitively affect the results of our analysis; it suffices 
that w±.o(xi ~ a) ~ 1 holds in some average sense. 

The formal solution of our initial condition problem 
now reads as 

w±(x,t) = f dx'K±(x,x,t)w±, (x), 

where the kernel K is defined through 

(±£ -dt + Difll) K ± (x, x, t) = 6(t)6(x - x). 

This equation is equivalent to Eq. (|i~4| ) for the propagator 
II and we can directly transcribe the solution. Substitut- 
ing 



K±{x -x, t) = T (27r)-^- iJ/2 (det a ± )-' /z x 



i(x-x(t)) CT± 1 (x-x(t)) 



with 



„±(A Q +A /3 )t _ 1 

a±(t) = ±2 ]T - ; ~ -v Q (v aT Dv> 3T 



a/3 



A Q + A/3 



into the integral representation for w± , one obtains 
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= T(2tt)- 



-(2d-l)/2 



(det a±) 



-1/2 



J dxe-^ x - s ( t » ff i 1 ( x - s W) U ; ± ,o(x) 



T(2tt) 



-(2d-l)/2 



(det ct±) 



-1/2 



/" dxe-^ x-i),7 ± 1(x -* ) «; ± ,o(x(-t)) 



With this background, we return to the discussion of the 
behavior of the generators w±(x!) on the reference tra- 
jectory. For times t < A -1 (corresponding to energies 
e > h\), <j± « u_f Aftl, implying that only points x in 
the imm edia te vicinity of the trajectory contribute to the 
integral (Bl). As a result, w±(x',t) <C 1 and the Green 
function remains essentially metallic. For t > A -1 , the 
component of the matrix tr + corresponding to the largest 
Liapunov exponent, behaves as (cr+)oo ~ vf^fs^ / 'A, 
and the integral begins to sample an exponentially grow- 
ing region around the reference trajectory. For 
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field configurations w+o(2/i ~ a ) = 0(1) contribute to 
the integral, which means that w+(x',t > £e) = 0{1). 
This identifies the inverse of the Ehrenfest time as the 
crossover scale beyond which the Eilenbcrger Green func- 
tion, even on exceptionally long trajectories, is locked to 
the superconductor order parameter. 
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